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Abstract 

We derive one-loop renormalization group equations (RGE’s) for Yukawa coupling pa¬ 
rameters of quarks and for the vacuum expectation values of the Higgs doublets in a general 
framework of the Standard Model with two Higgs doublets (2HDM “type III”). In the model, 
the neutral-Higgs-mediated flavor-changing neutral currents are allowed but are assumed to 
be reasonably suppressed at low energies. The popular “type II” and “type I” models are 
just special cases of this framework. We also present a numerical example for the RGE 
flow of Yukawa coupling parameters and masses of quarks. Detailed investigation, through 
RGE’s, of connection between the low energy and the high energy structure of these 2HDM’s 
would give us some additional insights into the viability of such frameworks, i.e., it would 
tell us for which choices of phenomenologically acceptable values of low energy parameters 
such frameworks can be regarded as reasonably natural (- without excessive changes in the 
parameter structure as energy of probes increases over a wide range). Furthermore, such 
analyses could provide us with possible signals of new physics at high energies. 

PACS: 11.10.Hi; 12.15.Ff; 12.15.Mm; 12.38.Bx; 12.60.Fr 


1 Introduction 

Low energy experiments show that flavor-changing neutral currents (FCNC’s) are very sup¬ 
pressed in nature. For example, 

Br(Kf — y + y~) — (7.4 ± 0.4) • 10~ 9 , — m B o | ~ (3.36 ± 0.40) • 10 -1 °MeV , 

\rriK L — niK s \ — (3.510 ± 0.018) • 10~ 12 MeV , \m D o — m D o| < 1.32 ■ 10 _1 °MeV , 

Br{b -> si) = (2.32 ± 0.67) • 10~ 4 etc. 
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The various alternative models of electroweak interactions - extensions of the minimal Stan¬ 
dard Model (MSM) - must take into account this FCNC suppression. The most conservative 
extensions of the MSM are apparently the models with two Higgs doublets (2HDM’s). The 
conditions for the one-loop FCNC suppression of contributions coming from gauge boson 
loops, i.e., the allowed representations of fermions, have been investigated some time ago ]l|. 
In addition to the MSM (one Higgs doublet model), Glashow and Weinberg Jl| proposed for 
the Higgs sector the “type I” and “type II” 2HDM’s. They proposed them as models which, 
in a “natural” way, have the zero value for the flavor-changing renormalized (i.e., low energy) 
Yukawa couplings in the neutral sector (called from now on: FCNC renormalized Yukawa 
couplings). These two types of the 2HDM’s have been widely discussed in the literature. 

Later on, extensions with more than one Higgs doublet, other than the 2HDM(I) and (II), 
have been proposed and investigated - the general “type III” 2HDM’s @-|2(J. This is the 
framework in which the renormalized (i.e., low energy) FCNC Yukawa coupling^ parameters 
are in general nonzero, but must be sufficiently suppressed. By “sufficiently” we mean that 
the measured low energy flavor-changing processes, including the A F = 2 processes, K°-K , 
B°-B and D°-D , do not get enhanced by tree level neutral Higgs exchange diagrams beyond 
the bounds given by low energy experiments. However, it may well be that the suppression of 
FCNC’s is not relevant for the top quark, i.e., we may have appreciable renormalized coupling 
of t and c to neutral scalars. This is so because, up to date, no stringent experimental bounds 
for these direct FCNC couplings of the heavy top quarks exist. 

The general 2HDM(III) framework has been first mentioned already in 1973 by T.D. Lee 
[|2|, for the hadronic sector and with the emphasis on the CP-violating phenomena originating 
from the nonzero relative phase of the two vacuum expectation values (VEV’s). Later on, 
the model has been investigated by several authors ||-|J, [jlTJ, who mainly investigated 
the bounds on the low energy (i.e., renormalized) parameters in the scalar mass and in the 
Yukawa parameter sector - the bounds resulting from available phenomenological data on 
CP-violating phenomena, such as e and s' parameters of the kaon physics, and the neutron 
electric dipole moment. 

In 2HSM(III), specific ansatze for the FCNC Yukawa coupling parameters^ have been 
proposed by Cheng, Slier and Yuan (CSY) 0, |£j, and by Antaramian, Hall and Rasin 
(AHR) 0. These two groups of authors also investigated bounds on the masses of scalars 
and on FCNC Yukawa coupling parameters arising from available phenomenological data 
of FCNC phenomena, such as the K-K and B-B mass differences and rare B decays. The 
basic messages of their works are the following: 

• their ansatze are reasonably natural (or, more cautiously: not “unnatural”) from the 
aspect of actual hierarchy of fermionic masses, since they are motivated to a large 
degree by this mass hierarchy. 

• their ansatze allow the masses of neutral scalars to be as low as rv./ 10 1 2 GeV while still 
not violating available (low energy) data of FCNC phenomena. 

CSY ansatz is given explicitly in Section 2. Later on, several authors further investigated 
implications of the CSY, AHR and/or related ansatze, and of specific assumed ranges of 

1 A more precise expression would be “neutral flavor-changing scalar (Yukawa) coupling,” since these 
couplings have no four-vector current structure involving 7 ^. 

2 Implicitly, all these coupling parameters are the renormalized (or: nearly renormalized) parameters, i.e., 
they are parameters at the low evolution energies since the authors investigated their tree level contributions 
to (low energy) FCNC phenomena for which phenomenological data are available. 
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masses of scalars, for the FCNC phenomena measured presently or to be measured in various 
possible future colliders, such as: (tree level effects in) the decays B —* /x + /x — 0; one-loop 


ey 


H u 


tc [|I 0 ; (t 


ree 


cZ |10|; two-loop effects in /i 
—> tc [|T7|; e + e~ —> Z* —> H°A° —> bbtc, W + W~tc, or ttcc |I3| ; 

mW- one-loop processes 77 —> h°, A 0 —> tc 


tc 18 


processes t —> 07 and t — 
level) process —> H° - 

one-loop process e + e~ —> 7 
[T4|1 ; e + e“ —> Z* — > (h°,H°)Z —> £cZ [Id]]; gauge boson (IUIU and ZZ) fusion processes 
e + e^ — > tcv e v e , tce + e~ |20fl ; rare decays t — > cLF+IU - , cZZ |20fl . 

Hall and Weinberg [U| pointed out that, while such 2HDM(I1I) models generically (nat¬ 
urally) suppress FCNC reaction rates to acceptable levels (e.g., by possessing approximate 
global flavor U( 1) symmetries leading to AHR-type of ansatze), they may in general give 
too much CP violation in the neutral kaon mass matrix. Therefore, they argued that these 
models must also possess CP as a good approximate symmetry. 

Low energy phenomenology of a special 2HDM(11I) framework, in which the renormalized 
FCNC Yukawa couplings (not involving the top quark) are suppressed practically to zero 
and the additional FCNC contributions come only from loops involving charged scalars, has 
also been investigated [^ 1 ] . 

In the present work, we construct RGE’s for the Yukawa coupling parameters of quarks 
(and for quark masses) in the discussed 2HDM(III) and thus obtain a means for investi¬ 
gations of the high energy behavior of such theories. The main motivation for the latter 
investigations can be summarized in the following question: For which (if any) phenomeno¬ 
logically acceptable ansatze of the FCNC Yukawa coupling parameters at low energies do 
we have a reasonable behavior of these parameters at higher energies of evolution? Under 
the “reasonable” behavior we understand a rather tame evolution of these parameters as 
the energy of probes increases by several orders of magnitude. Stated otherwise, this is the 
requirement that these parameters do not increase by order(s) of magnitude in the region 
of the evolution energy which is not very close (on the logarithmic scale) to the Landau 
pole]] of the top quark mass. Therefore, the theory would not change qualitatively in the 
FCNC Yukawa sector up until the energies of probes where the framework starts behav¬ 
ing nonperturbatively due to the large “mass” Yukawa coupling of the top quark. Such a 
reasonable behavior would then provide us with additional]] arguments that the discussed 
2HDM(III) frameworks, at least for some of the specific phenomenologically acceptable low 
energy choices of Yukawa parameters, are not unnatural. 

The Yukawa interactions in this 2HDM(III) framework in any SU(2) L -basis have the most 
general form 


- t {A?(sf * (1) )^ , + Dg , (5«*«)d« ) + 

i,j= 1 

+U-j\q ( 'l > + tj\f + /i.c.j + { UM-terms } . (1) 

The tildes above the Yukawa coupling parameters and above the quark fields means that 
these quantities are in an arbitrary SU(2) L -basis (not in the mass basis). The superscript (E) 
for the Lagrangian density means that the theory has a finite effective energy cutoff E, and 

3 Ref. Jl9| contains, in addition, an analysis of constraints from the electroweak p-parameter and from 
experimental data of various low energy processes such as A F = 2 processes, Z —> bb , rare B decays, 
e + e _ —> bs. 

4 Near the Landau pole of the top quark mass the theory starts behaving generally in a nonperturbative 
manner and the perturbative (one-loop) RGE’s start losing predictive power. 

5 - additional to the low energy arguments of CSY and AHR. 


r{E) 

^Y(III) ~ 
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the reference to this evolution energy E was omitted at the fields and at the Yukawa coupling 
parameters in order to have simpler notation (E ~ 10 2 GeV for renormalized quantities). 
The following notations are used: 


$( fc ) = 


/0dO+\ 

V 0 (*)° ) 


1 Yd + i#'1 
v/2 U l) + #i‘V 


= !r 2 4 (t| * 


1 / \ 


t(0 




($(!)) 




($ (2) ) 0 = 


e i?72 / 0 \ 

72 UJ ’ 


2,2 2 
V 1 + v 2 = V * 


( 2 ) 

( 3 ) 

( 4 ) 


In ([|), v [= v(E)] is the usual VEV needed for the electroweak symmetry breaking, he., 
v(E ew ) « 246 GeV. The phase difference g = r) 2 — between the two VEV’s in (|4j) may be 
nonzero; it represents CP violation originating from the purely scalar 2HD sector (cf. |L| ). 
The leptonic sector has been omitted in (jl[). 

We note that the popular “type I" and “type II” models are special cases (subsets) of 
this framework, with part of FCNC Yukawa coupling parameters being exactly zero 


2HDM(I): U {1) = L> (1) = 0 2HDM(II): V (1) = L> (2) = 0 . (5) 


In the 2HDM(II), the family symmetry in Cy enforcing this complete FCNC Yukawa param¬ 
eter suppression at all evolution energies E is of U(l)-type: d!'^ —> e l “<47 4>d) —> 
(j=l,2,3), the other fields remaining unchanged. This symmetry ensures that, in the course 
of renormalization, no loop-induced (In A cutoff-dependent) Yukawa couplings other than 
those of the form of the 2HDM(II) can appear. In the 2HDM(I), the family symmetry is 
similar: —> e ia d^\ u'^ —> e~ ia u^\ $^ 1 ' ) —>■ i n contrast to type I and type II, in 

type III 2HDM’s there is no exact (family) symmetry enforcing the complete suppression of 
the FCNC Yukawa coupling parameters. Stated otherwise, while FCNC Yukawa parameters 
in this general framework, at least those not involving the top quark, must be made quite 
small at low energies of probes E ~ E ew , they in general may increase when the energy of 
probes E increases. If they increase by order(s) of magnitude in the energy region well below 
the Landau pole, then such models should be regarded as rather unnatural - their behavior 
in the FCNC Yukawa sector is then drastically different at higher energies of probes (not 
just very near the Landau pole) from the behavior at low (electroweak) energies. 

In Section 2, we discuss relations between various notations for the scalar isodoublets and 
between various bases of the quark bleds, and the ensuing changes in the representation of the 
Yukawa coupling parameters - in order to highlight the suppression conditions imposed at low 
evolution energies on the FCNC Yukawa coupling parameters of the 2HDM(III). In Section 
3, we derive one-loop RGE’s of the Yukawa coupling parameters and of the scalar fields (and 
hence of their VEV’s) in the described 2HDM(III) model, with the purpose of investigating 
the behavior of the framework at evolving energies of probes E. In Section 4, we show and 
discuss one typical numerical example of the evolution of Yukawa coupling parameters in this 
general framework, in particular the evolution of the FCNC Yukawa parameters. The values 
of these parameters at low energies (E ~ E ew ) were chosen according to the CSY ansatz, 
and therefore they fulfill the FCNC suppression restrictions discussed in this Introduction 
and in Section 2. Section 5 is a summary of the results and conclusions. 
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2 Conditions of FCNC suppression at low evolution 
energies 

The Lagrangian density ([[])[] can be written in a form more convenient for consideration 
of FCNC Yukawa coupling parameters by redefining the scalar isodoublets in the following 
way: 


<F /(1) = (cos/3)$ (1) + (smP)e- iri <S>W , 
<f> ,(2) = —(sin/5)$ {1) + (cos/3)e - * ,?< F {2) , 


where 

V = r/2 - Vi 

and 

tan (3 = — => cos (3 — — , sin (3 — — 
V\ v v 

Therefore, the VEV’s of the redefined scalar isodoublets are 


e-i. (#'<»)„ = -L F) , <*«)„ = — f° 


V2\0; 


( 6 ) 

(7) 

( 8 ) 

(9) 


The isodoublet is therefore responsible for the masses of the quarks, and with its 
couplings to the quarks is responsible for the FCNC couplings, as will be seen below. The 
original Yukawa Lagrangian density (|l|) of 2HDM(III) can then be rewritten in terms of these 
redefined scalar fields as 


4m = - Z {4f > (vL*' ,1> )dR+<5<' ,) (?®i' (1, )flg > + h.c. } 
i,j= 1 

- Z {Ait*®*' 1 ”)® + « a (#«' (2) )+ h.c. } , (10) 

i,j=1 


where the Yukawa matrices G ^ and G^ are rescaled mass matrices, and U and D the 
corresponding “complementary” Yukawa matrices, in an (arbitrary) S'C(2)^-basis 


= = { C osf3)U {l) + {sin f3)e~ i71 U (2) , 

v 

G^ = —M(d) = (cqs/3)£) (1) + {sm/3)e +iv D {2) ; 
v 


( 11 ) 


U = ~(sm{ 3 )uW + (cos(3)e- iri uW , 

D = -{sin (3)D {1) + {cos f3)e +iri D {2) . (12) 

By a biunitary transformation involving unitary matrices , V® and V^, the Yukawa 

parameters can be expressed in the mass basis of the quarks, where the (rescaled) mass 

6 Throughout this Section we omit, for simpler notation, reference to the arbitrary evolution (cutoff) 
energy E at the quark fields, at the scalar fields and their VEV’s and at the Yukawa coupling parameters. 
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matrices and G^ are diagonal and real 


G<W = — M {u) 
v 

U 

G {D) = 

V 

D 


if G^lf * , M}P = S tj mP ; 

V?UVP ; 

VE& D '>V%' , Mg» = S,,mf ■ 
V?DVP | 


u L — Iff , u R — Vgu R , d L — Vgd L , — Vg d R . 


(13) 


(14) 

(15) 


The lack of tildes above the Yukawa coupling parameters and above the quark fields means 
that these quantities are in the quark mass basis (at a given evolution energy E ). The 
Lagrangian density (|T0|) can be written now in the quark mass basis. The neutral current 
part of the Lagrangian density in the quark mass basis is 


r (E) 

^■'Y(III) neutral 


-T E {Gif’dMW? 1 * + #?’) + 

+Gpu i l ) v$(<f>$ 1) - i<f>. 4 (1) ) + h.c. } 

—7S E {Dij^L^R^T +i(t>T ] ) + 

V 2 jj=i 


+U ij v$u%\</> 3 (2) - # 4 (2) ) + h.c. } . (16) 

On the other hand, the charged current part of the Lagrangian density in the quark mass 
basis is 


r {E) 

-‘-'Y(III) charged 


£ {(vG^) ij ufd'i\4 1) + ,4 1) ) - 

-(V'G^dgug^'P - #?>) + h.c. } 

-^£{ (VD^ufdPttP + id'P)- 


-(V U) tj dpu ( i\d>'} 2) - + h.c. } . 

Here, we denoted by V the Cabibbo-Kobayashi-Maskawa (CKM) matrix 

V = VfcKM = If V'" 1 . 


(17) 


(18) 


We see from (|IG| ) that the U and D matrices, as defined by (|I2D and (p~3|)-(pT^) through the 
original Yukawa matrices and D (j ' > of the 2HDM(III) Lagrangian density ([I]), allow the 
model to possess in general scalar-mediated FCNC’s. Namely, in the quark mass basis only 
the (rescaled) quark mass matrices G^ and G^ of (p~3D-(|T4|) [cf. also (pT|)1 are diagonal, 
but the matrices U and D in this general framework are in general not diagonal 


r(E) 

^Y(III) FCNC 


1 

2^2 


3 

E 

M =1 


{D + 03 2} +i(D- -D f )(d (i) d (i) ) 0 j 2) 
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+ (d + Z>t) 04 (2) - i (d - £>*).. (d^i^d^) 03 (2) 


1 

2^2 


£ [(" + "V 


u^u 


i,j=l 

*¥=3 



(u - C/ f ) (h ( V j) ) 



- (f/ + t/ f ) (uW^gu^) 04 (2) -i(u- U').. (u®i'y sU W) 03 (2) 


(19) 


It should be noted that the original four Yukawa matrices and (j = 1, 2) in an 
5'17(2)i-basis are already somewhat constrained by the requirement that (at low energy) the 
squares of the linear combinations^ and Mf°) are diagonalized by unitary transforma¬ 
tions involving such unitary matrices V^ and V ®, respectively, which are related to each 
other by = V. Here, V is the CKM matrix which is, for any specific chosen phase 

convention, more or less known at low energies. 

In order to have at low evolution energies (E ~ E ew ) a phenomenologically viable sup¬ 
pression of the scalar-mediated FCNC’s, the authors Cheng, Slier and Yuan (CSY) 0, 0 
basically argued that the elements of the U and D matrices (in the quark mass basis and at 
low evolution energies E) should have the form: 


Uij{E) 


c(u) V2 

j v 




Di,(E) = 


_ t (J)V2 


(rf) (d) 
m )j 


( 20 ) 


where 

eg°,d?~lfor E~E ew . (21) 

This form is in general phenomenologically acceptable and is motivated basically only by the 
actual mass hierarchies of quarks and the requirement that there is [at a given low energy of 
evolution (~ E ew )] no fine-tuning in which large Yukawa terms ufy (and: D^l) add together 
via flT2|) to make small terms Ujk (. Dj k ).| Therefore, this (CSY) form is considered to be 
reasonably natural. Similar (but not identical) ansatze have been proposed by the authors of 
[0] (AHR), motivated by their requirement that the Yukawa interactions have approximate 
1 /(1) flavor symmetries. 

From the CSY ansatz (|20|)-(f2T|) we see that the scalar-mediated FCNC vertices involving 
the heavy top quark are the only ones that are not strongly suppressed (at low evolution 
energies), since, as mentioned in the Introduction, FCNC processes involving the top quark 
vertices (not loops) are not constrained by present experiments. Later in Section 4 we will 
use low energy conditions for a numerical example of RGE flow of FCNC Yukawa 

coupling parameters. 


3 Renormalization Group Equations (RGE’s) in the 
general 2HDM 


7 Strictly speaking, the following “squares”: t and 

8 To visualize this point, Eqs. (Jllj) and ©) should be inspected, but this time in the quark mass basis 
(i.e., no tildes over the matrices). Then (0 would suggest that Uj l J {Dj 1 ^ ) is in general non-diagonal and 

of the order of \J/v (yj/v). As a result, ( |l2| ) would suggest that Ujk {Djk) is also of that 
order of magnitude, unless there is some peculiar fine-tuning on the right-hand side (RHS) of (0- For the 
complete suppression of FCNC Yukawa couplings ( Ujk = 0 = Djk for j ¥ k) we would need fine-tuning on 
the RHS of ©. 
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3.1 RGE’s for the scalar fields 


Here we present a derivation of the one-loop RGE’s for the scalar fields (||). The one-loop 
RGE’s for the Yukawa coupling matrices D^ k \ U^ (k = 1,2) will be presented in the next 
Subsection. In both derivations we will follow the finite-cutoff interpretation of RGE’s as 
discussed, for example, by Lepage . 

In order to calculate evolution of the scalar fields (p[ k \E) with “cutoff” energy E, we 
need to calculate first the truncated (one-loop) two-point Green functions —iT,[ k,t \p 2 -,E 2 ) 
represented diagrammatically in Fig. 1. More specifically, we calculate their cutoff-dependent 


(k) 


P 



(l) 



P 


Figure 1: The diagram leading to the two-point Green function E 2 ). Full lines represent quark 

propagators. 


parts oc p 2 In E 2 which will be ultimately responsible for the effective)] kinetic-energy-type 
terms ~ d u (p[ k \E)d u (p < f\E). In the course of the calculations, we ignore all the masses 
m ~ E ew of the relevant particles in the diagram. This would be consistent with the picture 
of a framework with a finite but large ultraviolet energy cutoff E E ew . For this reason, we 
don’t have to work in the mass basis of the relevant particles - these particles are regarded 
as effectively massless in the approximation, the transformations between the original bases 
of the relevant fields and their mass bases are unitary, and therefore the (mass-independent 
parts of the) calculated Green functions are the same in both bases. 

Calculation of the mentioned two-point Green functions — iE^’^(p 2 ; E 2 ), whose external 

(scalar) legs (p'^ and (jfp are truncated, is in the mentioned framework rather straightfor¬ 
ward. The relevant (massless) integrals over internal quark-loop momenta q can be carried 
out, for example, in Euclidean metric [q = (— iq 0 ,— g J ), p = (— ip 0 ,— p 3 ')], where the up¬ 
per bound in the loop integral is: q 2 < E 2 . After rotating back into Minkowski metric 
( p 2 i—> —p 2 ), we end up with the following results:^] 

1. Green functions with the external legs (p ( j k ' > and (pf^ having the same scalar indices 

(* = JY- 




V; £ 2 ) 


. K 




9 For clearer notation, we denote in this Section the evolving (UV cutoff) energy E at the fields not as a 
superscript, but rather as an argument. 

10 Throughout this Section, we implicitly ignore the cutoff-independent parts of the Green functions, since 
these parts are irrelevant for the mass-independent (M5-type of) RGE framework considered here. 
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where j = 1,2,3,4 (no running over the repeated indices j); k,£ = 1,2; m is an 
arbitrary but hxed mass of the order of electroweak scale (m ~ E e w ); and k stands for 


k = 


N c 

167T 2 


2. Green functions with the external legs (j>\ k ^ and (jE 1 having different scalar indices 






E 2 


41,2), 


K 


= ——n 2 In 
2 


Tr 


m- 


t/ (1) U (2) t _ u( 2 )(jW _ t (E) . (23) 


The above Green functions are antisymmetric under the exchange of the Higgs gener¬ 
ation indices k ^ £, and also under the exchange of the scalar indices i ^ j 


-ze^ 2) = +ix\y > = -^r = +*£& 

= -fS^ 2) = +fE^ 1} = -iY%f = 


42,1) _ 


42,1) _ 


(1,2) 


(24) 


3. The other Green functions are zero 


- = 0 


(25) 


for i = 1,2 and j = 3,4; for i — 3,4 and j = 1, 2; for i ^ j and k = i. 


All these Green functions can be induced at the tree level by kinetic energy terms. For exam¬ 
ple, in theory with the UV cutoff E , the kinetic energy term d u (j)[ k \E)d l/ (f)f\E) induces (at 

the tree level) the two-point Green function value (p 2 ; E 2 ) = ip 2 if (^(E) ^ <)Jp(E), 

and the value ‘lip 2 if (j)\ k \E ) = <p'p (E). Now, following the finite-cutoff interpretation of 
RGE’s as described, for example, by Lepage |23|, we compare the kinetic energy terms in 
the theory with the UV cutoff E and the equivalent theory with the slightly different cutoff 
(E + dE). The two-point Green functions in these two equivalent theories must be identical. 
When imposing this requirement in the tree + one-loop approximation, this leads to the 
following relation: 


1 4 

z t=i 
i 


4b 


bf 


1 0 
0 1 




3 =1 
K 


(E) 

(1) , , j>f]{E + dE) 


d v 


3 ’ ^3 

4 


1 o 
o 1 


(E) = 


d v 


bf 

i(2) 


[E + dE) 


(d\nE 2 )Y,d„ 

3 =1 


4i) 


1 (E) 


A n {E) A 12 (E) 
A 21 (E) A 22 {E) 


d v 


(dln£ 2 )£(-l)^ 


(hi) 


4>, 


<t>? 1 (E) 


0 B\ 2 {E) 

B 2 i(E) 0 


4 

4 2 ) 


d v 


( E) 

' 

A?) 


(E) , (26) 


where: the summation in the last sum runs over (i,j) = (1, 2), (2,1), (3,4), (4, 3); d(lnE 2 ) = 
In (E + dE) 2 — In E 2 = 2 dE/E\ and the real matrix elements A kf and B kk are related to the 
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one-loop two-point Green function expressions fl22|) and 


), respectively: 


Aki(E) — -Tr 


t _)_ u^)jj( k ) t _)_ jjWjjW t _j_ 


R 12 (F) = -Tr [G ( 1 ) G (2)t - G ( 2 ) G {1)t - Z^Z?^ + Z> ( 2 ) Z>^ 1 ) t j (F) = B 21 (E) 


(27) 

(28) 


Equation (g 6 |) is described in the following way: the sum on the left-hand side (LHS) and the 
first sum on the RHS represent the full kinetic energy terms of the scalars in the formulation 
with the UV cutoff E and (E + dE), respectively. The one-loop contributions of Fig. 1 with 
the loop momentum \q\ in the energy interval E < \q\ < A are already contained in the kinetic 
energy terms of the LHS effectively at the tree level (A is a large cutoff where the theory 
is presumed to break down). On the other hand, the kinetic energy terms of the (E + dE) 
cutoff formulation [the first sum on the RHS of (|26j)] effectively contain, at the tree level, 
the one-loop effects of Fig. 1 for the slightly smaller energy interval: (E + dE) < \q\ < A. 
Therefore, the one-loop two-point Green function contributions^ —idT,[ k j £ \p 2 -,E 2 ) of Fig. 1 
from the loop-momentum interval E < \q\ < (E + dE) had to be included on the RHS of 


(|26|) - these are the terms in the last two sums there. This is illustrated in Fig. 2. 


4 / 

jdq 


i,j,k,l 


(271) 

E<lql<E+dE 



+ 


(l) 


,(k) , , ,(l) 

^ 4>j 


= 0 


Figure 2: Diagrammatic illustration of the RGE relation (|2^) leading to the evolution of the scalar fields. 
<pZ stands for cj)[ k \E), and dcf> stands for (f>(E + dE) — cf>(E) (</> is a generic notation for The cross 

represents the contribution of the change of the kinetic energy terms originating from the changes d(f> of 
scalar fields. 

In order to find RGE’s for the scalar fields 4>[ k \E), we make the following ansatz for a 
solution of Eq. 


4>i{E + dE) < = 


t(2) 


(E + dE) > = <j>i{E) + da(i- E)4»i{E) + d/3(r, E)Z(E) , (29) 


where da{i\ E) and d/3(i] E) are infinitesimally small 2x2 matrices 

doty i da 


da(i; E) = 


Z Z 


0 d/3 12 

dp 2i 0 




(30) 


and the index i! in (^9|) is complementary to index i: 

i = 1^ i' = 2 , i — 2 i-> i! — 1 ; i — 3 i' — 4 , i — 4 i-> i' — 3 . 

Inserting ansatz (]29]) into the RGE relation (^6]) , we end up with the following set of relations: 

da ke (i] E) + da e k{i; E) = -n (dlnE 2 ^ A ke (E) , (31) 

11 More precisely: the corresponding effective kinetic energy terms. 
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W l2 { 1; E) + d0 21 (2; E)\ = - [d/3 21 ( 1; £) + d0 12 (2; £)] = -k (din E 2 ) B 12 (E) = 

= [d(3 12 (3-E) + dP 21 (4-E)] = ~[df3 21 (3-E) + d(3 12 (A-E)] . (32) 

In principle, these relations alone do not define the elements da k ^{i\ E) and d(3 ki (i] E). How¬ 
ever, there is another requirement that should be imposed on these transformation coeffi¬ 
cients: the resulting RGE evolution of the isodoublet fields and & 2 \E) should be 

invariant under the exchange of Higgs generation indices 1 2, because these two Higgs 

doublets appear in the original Lagrangian density ([I]) in a completely 1 2 symmetric 

manner. We will see in retrospect that this discrete symmetry is respected once we impose 
the conditions 



da M (i] E) 

= da lk {i\ E) , 

(33) 

df3 12 (l;E) 

= df3 21 (2,E) 

, d(d 12 (2; E) = df3 21 (l, E) , 


d(3 i 2 (3; E) 

= d/3 21 (4,E) 

, d(3 12 (4;E) = d(3 21 (3,E) . 

(34) 


Solutions (M)-(M) of the RGE condition (|26|), together with the symmetry conditions 
(|34]), lead to specific expressions for the evolution coefficients da.ki(i] E) and d/3 ke (i ; E). When 
inserting these coefficients back into the scalar field evolution ansatz (|29|), we end up with 
the following RGE’s for the evolution of the scalar fields: 


167T 2 d 
N c dlnE 

—-Tr 
2 


$\E) = -Tr fc/Wc/Wf + 0^ 


t/(i)C/(2)t + 2 )t + D (2) D (1) t] 0 


( 2 ) 


+i(-l) J -Tr 


[yf 1 )[/( 2 )t — E/f 2 )t/Glt _ _|_ jjWjji 1 )t 


02 ) 


(35) 


167T 2 d 


(^f\E) = -Tr [H (2) H (2)t + Z) (2) Z) (2)t l 0 


N r dlnE 2 


( 2 ) 


-Tr 

2 


C/(i) ff( 2 )t + H (2) H (1)t + Z) (1) l) {2)t + D (2) D (1)t 1 0j 1} 


+i(-iy +1 -Tr 


f/(i)C/(2)t - H (2) H (1)t _ £>(l)£)(2)t + ^(2)^(1)! 0(1) ) (36) 


where again / is the scalar index complementary to index j\ ( j,j') = (1,2), (2,1), (3,4), 
(4,3). These RGE’s lead to RGE’s for scalar isodoublets 


16n 2 d 
N c dlnE 


4 > (fc) (£) 


qy jjl k )jjd)l _)_ i)d)jj{ k )l <j>h) 


r=i 


(37) 


We see indeed that this set of one-loop RGE’s is invariant under the exchange 1 <->• 2, as 
required by the form of the Yukawa Lagrangian density (Jl|) of the 2HDM(III). 

In addition to quark loops, there are also loops of the electroweak gauge bosons contribut¬ 
ing to one-loop two-point Green functions of the scalars. However, since these gauge bosons 
couple to the Higgs isodoublets identically as in the minimal Standard Model (MSM), their 
contributions to the RHS of RGE’s (|35| ) - (|36|) and ([77]) are the same as in the MSM.j^] Con¬ 
sequently, the full one-loop RGE’s for the evolution of the scalar isodoublets in 2HDM(III) 

12 For these contributions of EW gauge bosons in the MSM, see for example Arason et al. |24j], App. A. 
However, note that they use for the U{l)y gauge coupling g\ a different, GUT-motivated, convention: 

(ffl)Arasonet al. = (5/3)(ffi)here- 
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are 


167r 


d 


d In E 


4> ( k \E) = -N C J2 Tr [c7 (fe) C/ (£)t + D^D (k) f ] 4> w 
£=1 

3 2 / T 71 \ , 9 2/ 


^ 1 2 (£) + ^ 2 2 (£) 


4 > (fc) (£). 


( 38 ) 


These RGE’s are simultaneously also RGE’s for the corresponding VEV’s 


167 T 2 —p— (e ir,k v k ) = —N c Tr \u {k) U w + D^D^] 


d In E 


+ 


i =i 
2 / 


Iflft-E) + jpKD 




(39) 


In this paper we don’t discuss the question of quadratic cutoff terms A 2 which appear in 
the radiative corrections to VEV’s in any SM framework. In the MSM, their consideration - 
under the assumption of the top quark dominance of the radiative corrections in the scalar 
sector - leads to severe upper bounds on the ultraviolet cutoff A for a substantial subset of 
values of the bare doublet mass and of the bare scalar self-interaction parameters M 2 (A) 
and A(A) - cf. Ref. [p5f]. 


3.2 RGE’s for the Yukawa matrices 

In order to derive one-loop RGE’s for the Yukawa matrices U^ and D^ k \ we will need the 
results of the previous Subsection concerning evolution of the scalar fields. In addition, we 
will need evolution of the quark fields u and d^ R - The latter can be derived in close 
analogy with the derivation of the evolution of scalar fields of the previous Subsection. Now, 
the diagrams (Green functions) of Figs. 1 and 2 are replaced by those of Figs. 3 and 4, and 
the scalar field kinetic energy terms in ( |2B| ) are replaced by those of the quark fields. The 
one-loop two-point Green function of Fig. 3, with the incoming and outgoing u ( ' ] '> of 


Ji) 

u 



Jj) 

u 


p p-q p 

Figure 3: The diagram leading to the two-point Green function — iE(p; E\ u(*),uW)). Dashed and full lines 
represent scalar and quark propagators, respectively. 


momentum p, in the framework with UV cutoff E, is 


-iS (p;E-,u®,u&) 


in (I?) 4 (i - 7s) £ 

+2(1 + 75 )]T . 


( 40 ) 
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The Green function with the incoming and outgoing of momentum p is obtained 
from the above expression by simply exchanging U^' 1 and . We now 

make the ansatz for the running of the quark fields 

duW(E) LjR { = u {k \E + dE) LiR - u (k \E) L ^ = df u (E)i L e ’ R) uW(E) L , R , (41) 

dS k \E ) LtR | = d (k \E + dE) L)R - & k \E ) LtR | = df d (E)i L d R) d^(E) L , R , (42) 

where subscripts L, R denote handedness of the quark fields: q L = (1 — 75 ) 9 / 2 , q R = 

(1 + 7 5 )g /2 (q = u^ k \dfE>). In complete analogy with the previous Subsection, we obtain 

from these ansatze and from the RGE relation illustrated in Fig. €3 the following relations 



Figure 4: Diagrammatic illustration of the RGE relation leading to the evolution of quark fields. Physically, 
this relation means that the two-point Green functions with truncated external (quark) legs, at one-loop 
level, are the same in the theory with E cutoff and in the theory with E + dE cutoff. Conventions are the 
same as in previous figures. 


for the quark field evolution matrices df u : 


«.(£)«’* + dk(E)f 

dfu(E)\f >* + df u (E)f 


32tt 2 


32t r 2 


(dlnE 2 ) Y \fjWuW + (E) , 

L J ji 


k=1 
2 


(dlnE 2 )y [u m il {k) ].. (E) . 


k=1 


(43) 

(44) 


The relations for the df d evolution matrices of the down-type sector are obtained from the 
above by simple exchanges and . A solution to all these relations 

for the quark held evolution matrices is 


dUE)f = 


y [umw 

k=1 

+mmP.xE) = df d (E)f , 

j ij j 

(45) 

dU(E )!f = 

-64 U dl » E2 ' 

y [u^u^y 

k=1 L 

..(E) , 
y 

(46) 

*(B)Sf = 

~y^ E2 ' 

y 

k=1 

( E ) • 

y 

(47) 


Another Green function needed for the derivation of the one-loop RGE’s of Yukawa 
matrices is the one represented by the diagram of Fig. 5. When the external legs there are 
(incoming, with momentum k ), u^) (outgoing, with momentum p + k ), and dp [or dP], 

13 The RGE relation represented by Fig. 4 is analogous to relation ( p6[ ) represented diagrammatically in 
Fig. 2, but this time the kinetic energy terms are those of the quark fields. 
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k 


Figure 5: One-particle-irreducible (1PI) diagram contributing to the evolution of the Yukawa coupling 
parameters. Conventions are the same is in previous figures. 

it turns out that only the diagram with the charged scalar exchange contributes, and the 
resulting truncated three-point Green function, in the framework with the UV cutoff E , is 

G <3) (*,p; 


The corresponding Green function with the down-type quark external legs is obtained from 
the above expression by the simple exchanges U ^ and . 

Now, the one-loop RGE’s for the Yukawa matrices are obtained in analogy with the rea¬ 
soning leading, in the case of two-point scalar Green functions, to the RGE relation (^) in 
the previous Subsection [cf. Fig. 2], It is straightforward to check that the contribution of 
the quark loops in the scalar external leg cancel the contributions coming from the renormal¬ 
izations of the scalar fields in the kinetic energy terms of the scalars - this is illustrated in 
Fig. 6. Furthermore, it can be checked that the contributions of the scalar exchanges on the 





327t 2 v / 2 


In 


E 2 


m- 


(1 + 7s) 


f)(r) 


Ajl 


+ (1 - 7s) 


jj(r) t £)(7 jj(r) t 




(48) 



E<lql<E+dE 


Figure 6 : Cancelation of contributions from the quark loop (one-particle-reducible - 1PR) with those of 
the scalar field renormalizations in the kinetic energy term of the scalars, for the energy cutoff interval 
(. E,E + dE). 

external quark legs cancel the contributions coming from the renormalizations of the quark 
fields in the kinetic energy terms of the quarks - this is illustrated in Fig. 7. All in all, the 
1PR one-loop contributions are canceled by the contributions of held renormalizations in 
the kinetic energy terms. Therefore, the only one-loop terms contributing to the evolution 
of the Yukawa matrices are those depicted in Fig. 8. The three diagrams with crosses 
there correspond to contributions of the following changes in the Yukawa coupling terms: 

• Yukawa matrix change (renormalization) dU^ [= U^ k \E + dE) — U( k \E)\ - Fig. 8(b); 
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0 



E<lql<E+dE 


Figure 7: Cancelation of contributions from the scalar exchange on the quark legs (1PR) with those of 
the quark field renormalizations in the kinetic energy term of the quarks, for the energy cutoff interval 
(E,E + dE). 



E<lql<E+dE 

Figure 8 : Diagrammatic representation of the RGE for the up-type Yukawa matrix U. Only the 1PI scalar 
exchange [(e)] and the effects of the renormalizations of the Yukawa matrix, of the scalar fields and the 
quark fields in the Yukawa couplings [(b), (c), (d), respectively] contribute when the cutoff is changed from 
E (RHS) to E + dE (LHS). Note that dU stands for U(E + dE) — U(E), etc. The contributions of the gauge 
boson exchanges were not considered in the Figure. 


• the scalar field renormalization dcf)^ [= 4>^\E + dE ) — 4>^(E )] - Fig. 8 (c); 

• the quark field renormalization du^ and du^> - Fig. 8 (d). 


Figure 8 is a diagrammatical representation of the physical requirement that the three-point 
(quark-antiquark-scalar) Green function, at one-loop level, be in the theory with the cutoff 
E + dE [left-hand side of Fig. 8 : (a) + ... + (e)] the same as it is in the theory with the 
slightly lower cutoff E (right-hand side). 

Using the results of this and the previous Subsection, we can then write down the one- 
loop RGE for U^ corresponding to Fig. 8 , at the right-handed component [oc (1 + 75 )] of 
the three-point Green function 


+ dU$ + ^— 2 (dhiE 2 ) I -iY c ]TTr [u^U™ + U® 


1 = 1 


— £ t) U™ + 2 uWuMU® 


1= 1 


+2J2 


= u {k) 


(49) 


t=\ 




The first sum on the LffS (oc N c ) corresponds to Fig. 8 (c) [cf. Eq. (|37|)1, the second sum to 
Fig. 8 (d) [cf. Eqs. (ff5p, (|46|)1, and the third sum to Fig. 8 (e) [cf. Eq. (^BD). The left-handed 
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part of the Green function yields just the Hcrmitean conjugate of the above matrix relation. 
The analogous consideration of the three-point Green functions with the down-type external 
quark legs and dP^ gives relations which can be obtained from the above relation again 
by the simple exchanges <->• D l ' s ' > and . These relations can be rewritten in 

a more conventional form 


d 


16vr 2 ^— U {k \E) = <N c y Tr U ik) U {e ^ + D (e) D m U w 
d\nE v ’ ' ^ 


t =1 


.1 J2 'u®u ®t + D w £)W U {k) + U {k) y U {i) ^U (i) 


i =i 

-2y \pd)E>mij{t) 

e=i 


i =i 


(50) 


and an analogous RGE for D^ k \ These RGE’s still don’t contain one-loop effects of exchanges 
of gauge bosons. However, since the couplings of quarks and the Higgs doublets to the gauge 
bosons are identical to those in the usual MSM, 2HDM(I) and 2HDM(II), their contributions 
on the RHS of the above RGE’s are identical to those in these theories. Therefore, the final 
form of the one-loop RGE’s for the Yukawa matrices in the general 2HDM(III) now reads 


d 


16 7 r 2 -^—^G (fc) (E) = \ N c ~y Tr [c/ (fe) C/ wt + L> w L> (fe)t [ U® 


e= i 


.1"f/ w G wt + D (i) U {k) + U {k) y 


i =i 

-2y 

1 =1 


t =1 


- A v m , 


(51) 


d 


16n 2 —-— D^ k \E) = { N c Y Tr \D (k) D + U w U {k) H D w 
d In E 11 1 


i =i 

1 2 


'uWfjw + f)(k) + f)( k )y 


t=\ 

2 


*=1 


- 2 ^ 


£=i 


- A d D^ 


(52) 


where the functions H[/ and characterizing the contributions of the gauge boson ex¬ 
changes, are the same as in the MSM, 2HDM(I) and 2HDM(II) 


, 0 (iV c 2 -l) 2 9 2 17 2 

Atj — 3 - CJo + —(7o -f- Q | , 

7V C 3 4 2 12 1 

— g 2 , 

and the gauge coupling parameters g 3 satisfy the one-loop RGE’s 

d 


16vH 


d In E 


9j Pj9j ’ 


(53) 


(54) 
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with the coefficients Cj being those for the 2HDM’s (iV# = 2) 


C 3 = -(llN c -2n q ) , C -2 — 7 — -n q . G 1 = --- 


10 

Y 


-n„ 


(55) 


Here, n q is the number of effective quark flavors - e.g., for E > m t we have n q ~ 6; for 
77ib < E <m t we have n q ~ 5, etc. 

The obtained relevant set of RGE’s for the VEV’s e irik Vk ( |39| ) and the Yukawa matrices 
UW and (|Y|) -(|52D determines in principle also the running of the quark masses. Instead, 
we can rewrite all these RGE’s in a representation involving the VEV parameters v = 
Jvi + v%, tan (3 = v 2 /v\ and rj = r] 2 — rft [cf. (||)], and matrices & u \ & D \ U and G [cf. (|lT|), 
(HH)] ~ this representation is more convenient for discerning the running of the quark masses 
and of the FCNC couplings. Applying lengthy, but straightforward, algebra to the hitherto 
obtained RGE’s results in the RGE’s of the latter set of parameters 


I 67 G 


d 


d In E 




= -21Y.Tr 


q{U)q{U)\ + G(D)Qm] v * + 


3 2 9 2 
2 dl + 2 ^ 2 . 


(56) 


167T' 


d 


d In E 


(tan/3) = 


N r 


2 cos 2 f3 


Tr 


UG m + G (D) G f 1 + Tr \& u) U ] + GG (D)t 1 1,(57) 


1 a 2 d ( \ N, 
16vr — —( 77 ) = — 


dlnE 


i sin(2/3) 


Tr 


G {u) U ] - UG m ' 


Tr 


G (d) D f - DG [D)] 1 l , (58) 


16tA 


d 


dlnE 
1 


(U) = N c \ 2Tr \lJG m + & d) D'] G (u) + Tr \UU f + GG+ 


U 


(cot /3) Tr 


U + g (c°t /5)Tr 


G (u) U f - UG m 

UU t + D& + GMGW + G^G^ 
-2DD ] U - 2 G (d) D^G {U) - AuU j , 


G^G* - DG m 


U 




u + u \uW + GW'GW 


(59) 


167T 


d 


dlnE 


(G) = wJ 2Tr \DG m + G (U) U f l G (D) + Tr \lJU ] + TT 1 


D 


--(cot^Tr 


D + -(cot YTr 


G (D) G f - GG (S)t 

GC/t + DD' + gWgW + G^G™ 
2 UU ] D - 2G (c/) G t G (D) - A D G j , 


G^f/t - U&M f 


G 




G + D \D^D + G (D)t G (D) 


(60) 
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d 


16tt 2 —-- (G (r/) ) = NA Tr [G (c/) G (c/)t + & D) G m ] G (c/) 


dln£ 

-^(tan/3)Tr [G^t/ 1 - UG m ] G {U) + ^(tan/3)Tr [G {d) D ] - £>G (D)t ] G ([/) 


i[ 


UW + DD t + G^G^t + G^G^tl ^ + GW \fflU + G^G^ 


-2DG m U - 2G (D) G (D)t G ([/) - ^G^j , 


( 61 ) 


d 


1671 "^ ———- (G^) = 7V C Tr [G^G^t + G^G™] G (D) 


d In E 

-i(tan/5)Tr [g ( d) ^'- DG (D)t ] G (D) + ^(tan/?)Tr [G^C/t - f/G ([/)t ] G (D) 


2 L 


f/f/ f + ZXD 1 ' + G {u) G m + G (D) G (D)t l G (d) + G {d) \D ] D + G (D)t G (D) 


-2 UG m D - 2G (l7) G (c/)t G (D) - A D G {D) i . 


(62) 


4 A numerical example of evolution 

Here we present one simple but hopefully typical example of the RGE evolution of parameters 
in the 2HDM(III) framework. For simplicity, we assumed: 

• there is no CP violation - all original four Yukawa matrices U^\ are real, and the 
VEV phase difference 7] is zero; 

• the masses and Yukawa parameters of the first quark generation as well as that of the 
lcptonic sector are neglected (the quark Yukawa mass matrices are therefore 2x2). 

For the boundary conditions to the RGE’s, at the evolution energy E = Mz , we took 
the CSY ansatz (|20| ) -(|2l| ) , with = 1 = for all i,j = 1,2 [note: i — 1 refers now 
to the second quark family (c,s), and % — 2 to the third family (t,b)]. For the (2 x 2) 
orthogonal CKM mixing matrix V we took V 12 (Mz) = 0.045 = — V 21 {Mz)- The values of 
other parameters at E — Mz were chosen: 
tan/? = 1.0; v = \Jv\ + uf = 246.22 GeV; 
a 3 = 0.118, a 2 = 0.332, a x = 0.101; 

m c = 0.77 GeV, m s = 0.11 GeV, m b = 3.2 GeV, and m t = 171.5 GeV. 

The latter values of quark masses correspond to: m c (rri c ) ~ 1.3 GeV, m s (lGeV) ~ 0.2 GeV, 
m b (m b ) ~ 4.3 GeV, and m b hys ' ~ 174 GeV [m t (m t ) ~ 166 GeV]. For a :i (E) we used two- 
loop evolution formulas, with threshold effect at f? ~ m b hys ' taken into account; for otj(E) 
(j = 1 , 2 ) we used one-loop evolution formulas. 

The described simplified framework resulted in 18 coupled RGE’s [for 18 real parameters: 
v 2 , tan (3, Uij , b l]: G^\ Gyf'*], with the mentioned boundary conditions at E — Mz- 
The system of RGE’s was solved numerically, using Runge-Kutta subroutines with adaptive 
stepsize control (given in J2(J). The numerical results were cross-checked in several ways, 
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including the following: FORTRAN programs for the RGE evolution and for the biunitary 
transformations were constructed independently by two of the authors (S.S.H. and G.C.), 
and they yielded identical numerical results presented in this Section. 

The results for the FCNC Yukawa parameter ratios Uij(E) {Uij(M z ) and Dij(E)/Dij(M z ) 
(i 7 - j) are depicted in Fig. 9. From this Figure we immediately notice that the FCNC 



Figure 9: Evolution of the FCNC Yukawa parameter ratios Uij(E)/Uij(M z ), Dij(E)/Dij(M z ) (i A j) in 
2HDM(III). These parameters are in the quark mass basis. The choice of parameters of the model at the 
starting low energy E = M z is specified in Sec. 4. Yukawa couplings of the first generation were neglected; 
*=1,2 correspond to the second and third quark generation, respectively. 

coupling parameters of the down-type (b-c) sector are remarkably stable as the evolution 
energy increases. Even the up-type FCNC ratios, although involving the heavy top quark, 
remain rather stable. Only very close to the top-quark-dominated Landau pole (E po i e ~ 
0.84 • 10 13 GeV)[3 the coupling parameters start to increase substantially. For example, 
in the down-type FCNC sector (b-c) the corresponding ratio D 21 (E)/D 21 (M Z ) acquires its 
double initial value (i.e., value 2) at E & 0.7-E po i e , which is very near the (Landau) pole.Q 
For the ratio Di 2 (E) / D V2 (M Z ) the corresponding energy is even closer to E pole . 

One may ask whether the mentioned stability features in the evolution of all Yukawa 
coupling parameters in the framework. This doesn’t seem to be the case. For example, 
in Figs. 10 and 11 we depicted, with the described case of initial conditions, evolution of 
the Yukawa coupling ratios connected with no FCNC’s: Ujj(E)/Ujj(M z ), Djj(E) / Djj(M z ), 
G ( fj\E)/G < jp(M z ) and G^\e)/G^\m z ) (j = 1,2). We see that evolution behavior of 
many of these ratios, in stark contrast to the FCNC parameter case of Fig. 9, is far from sta¬ 
ble. For example, the ratio D U (E)/Du(M z ), connected with s quark, acquires double its ini- 

14 The value of F po i e is strongly dependent on the given value of parameter £, as shown later in Fig. 13. 

15 Approximately the same is true also for U 2 i{E)/U 2 \{M z ). 
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Figure 10: Same as in Fig. 9, but for the neutral current Yukawa coupling parameters Ujj and Djj (j = 1 , 2) 
which don’t change flavor. 

tial value approximately at the value 0.1E po i e ; the “mass” Yukawa ratio G < ^\e)/G < ^\Mz)- > 
corresponding to b quark, acquires half its initial value at 10~ 5 E po i e . These energies are 
therefore substantially further away from the Landau pole than the corresponding energies 
for the FCNC Yukawa down-type sector (cf. previous paragraph), even on the logarithmic 
scale. 

In Fig. 12 we depicted evolution of the quark masses for the discussed case. Note that 
they are not simply proportional to the “mass” Yukawa parameters G < fp{E) and G^\e), 
because the VEV v(E) also evolves with energy [cf. Eq. (|56|)]. Logarithmic scale was chosen 
for the masses in order to include rrit(E) in the figure. From this figure we see that we 
have a rather strong variation of rrib{E) [and also of mt(E)\ continuously when the evolution 
energy increases. For example, rrib{E) reaches half its initial value [rrib(E) = rrib(Mz)/2] at 
E ps 10 _6 E' po i e , which is very far away from E pole - this again contrasts with the behavior of 
FCNC Yukawa parameters of Fig. 9. 

The discussed numerical example of the 2HDM(III) framework tells us that there def¬ 
initely exist choices of reasonably suppressed (he., phenomenologically acceptable) FCNC 
Yukawa coupling parameters at low energies such that these parameters remain largely un¬ 
changed (suppressed) up to energy regions very close to the Landau pole. On the other 
hand, this behavior doesn’t feature in the entire sector of the Yukawa coupling parameters. 

It should be stressed that these results are independent of the chosen value of the VEV 
ratio tan/3 at E — Mz■ This is connected with our choice of the CSY boundary conditions 
0-0 at E — Mz for the Yukawa matrices in the quark mass basis (^“’ = = 1) 

and the reality of the chosen CKM matrix at E — Mz- These boundary conditions result 
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in real and /3-independent Yukawa matrices U, D, & u \ G ^ in a weak [SU(2)l] basis^ at 
E = M z . The RGE’s (|5P|)-(|^|) then imply that these matrices remain real and independent 
of /3 at any evolution energy E, and that also their counterparts U, D , and G^ in 
the quark mass basis, as well as the CKM matrix V, remain real and independent of (3 at 
any energy E. Stated otherwise, if there is /3-independence and no CP violation (neither in 



Figure 11: Same as in Fig. 9, but for the “mass” Yukawa parameters G^p and G^ (j = 1 , 2 ) instead. 
Since G^ U \E) and G^ d \E) matrices are diagonal by definition (quark mass basis), these neutral current 
Yukawa coupling parameters have zero FCNC components automatically. 


original Yukawa matrices nor in the scalar sector) at a low energy (E = M z ), then these 
properties persist at all higher energies of evolution.^ 

This feature is in stark contrast with the situation in the 2HDM(II) where the Yukawa 
matrices strongly depend on /3 already at low energies - e.g., g t (M z ) = m t (M z )\/2/v u = 
m t (M z )V2/[v sin(/3(M z ))\. Also the location of the Landau pole in the 2HDM(II) then 
crucially depends on (3(M Z ) - smaller /3(M Z ) implies larger g t (M z ) and hence a drastically 
lower Landau pole. 

On the other hand, the 2HDM(III) framework treats the up-type and the down-type 
sectors of quarks (the two VEV’s V\ and V 2 ) non-discriminatorily. Therefore, it should be 

16 We chose at E = M z the following weak basis: U = U, & u ' 1 = G^ u \ D = VD, = VG^ D \ where 

V is the CKM matrix (at E = Mz)- According to relations && the reality of the Yukawa matrices U, 
D, anc l G( d '> at the low energy E = M z would follow, for example, from: the requirement of no CP 

violation in the Yukawa sector (i.e., the original Yukawa matrices and ZF-b are all real) together with 
the requirement of no CP violation in the scalar sector (i.e., the VEV phase difference 77 = 0) at that low 
energy. 

17 CP conservation in the pure scalar sector at a low energy E = M z (i.e., 77 = 0) also persists then at all 
higher energies of evolution, since dr]/dinE = 0 by the reality of the Yukawa matrices, according to RGE 

©■ 
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Figure 12: Evolution of the quark masses m^ l \E) = G^-\E)v(E) / \/2 and (E) = G^\E)v(E)/V2 


( (d) 

{m\ = m s 


Ad) 


= mb ; rrif 1 = m c , = mi), for the discussed numerical example. 


A u ) 


expected that any reasonable boundary conditions for Yukawa coupling parameters at low 
energies should also be independent of /3 in such frameworks, and this independence then 
persists to a large degree also at higher energies. Also the locations of the Landau poles (i.e., 
of the approximate scales of the onset of new physics) should then be expected to be largely 
/1-independent. In this sense, 2HDM(III) has more similarity to the minimal SM (MSM) 
than to the 2HDM(II). The persistence of complete /^-independence of the Yukawa coupling 
parameters at high energies and of the Landau poles, however, can then be “perturbed” 
by CP violation - because RGE’s (|59|)- (02|) are somewhat /3-dependent when the Yukawa 
matrices U, etc., are not real. Also the VEV phase difference r) is then not a constant when 
the energy increases [cf. (|58|)1. 

In addition to the connection between (low energy) CP violation and /3-dependence of 
high energy results, there is yet another feature that distinguishes the 2HDM(III) framework 
from the MSM - the Landau pole of a 2HDM(III) framework is in general much lower than 
that of the MSM. We can see that in the following way: let us consider that only the 
Yukawa parameters connected with the top quark degree of freedom are substantial, i.e., 
G 2 P = g t ~ 1 and U 22 = g\ ~ 1. We have: g t (E ) = m t (E)\/2/v(E), as in the MSM, 
and g' t (E) is an additional large Yukawa parameter - both crucially influence location of the 
Landau pole. Inspecting RGE’s (|5U|) and ((JI) for this special approximation of two variables 
g t and g[, we see that RGE for g t is similar to that in the MSM, but with an additional 
large positive term on the right: (3/2 ){g' t ) 2 g t . RGE for g[ has a similar structure as RGE for 
g t , but with substantially larger coefficients at the positive terms on the right. As a result, 
g[[E) is in general larger than gt(E). Our specific numerical example [cf. Figs. 10 and 11 
for U 22 and G 2 P] shows that g[{E) is on average (average over the whole evolution energy 
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range) almost twice as large as gt(E). If we then simply replace in the mentioned additional 
term (3/2 )(g' t ) 2 gt the parameter ( g[) 2 by 3.5 g%, we obtain from the resulting “modified” 
MSM RGE for g t a value for the Landau pole in the region of 10 12 — 10 13 GeV, which is 
roughly in agreement with the actual value of the Landau pole of our numerical example 
Ep 0 i e ~ 0.84 • 10 13 GeV. And this value is much lower than F po i e in the MSM which is above 



Figure 13: Variation of the Landau pole energy when the low energy (E = Mz) parameters = £ 

of the CSY ansatz (|2^)-(|2^) are varied. For £ = 2.5, the onset scale of new physics is already quite low: 
E pole ~ 2 TeV. 


the Planck scale. Of course, when we allow the parameters of the CSY ansatz (|20|) - 
( PI| ) at E = M z to deviate from 1, we obtain larger log(F , pole ) for smaller and smaller 

log(F po ie) for larger Qj. In Fig. 13 we depicted this variation of the Landau pole energy 
when the CSY low energy parameters are varied. 


5 Summary and Conclusions 

We have derived the one-loop RGE’s for the quark Yukawa coupling matrices and the VEV’s 
in the Standard Model framework with the most general two-Higgs-doublet Yukawa sector 
[2HDM(III)]. A simple - and at low energies phenomenologically acceptable - numerical 
example for the resulting evolution of these Yukawa parameters suggests that the framework 
cannot be dismissed as “unnatural” from the FCNC-RGE point-of-view. Stated otherwise, 
the numerical example shows remarkable stability of the suppressed FCNC Yukawa coupling 
parameters when the energy of probes increases continuously all the way to the vicinity of 
the top quark Landau pole. The Landau pole is in general well below the Planck scale in this 
framework. We believe that further numerical investigations are warranted, in order to see 
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whether and/or to what degree this behavior survives when we scan over certain reasonable 
(phenomenologically acceptable) ranges the values of low energy parameters of the model - 
i.e., the Yukawa coupling parameters and tan (3 at E ~ E ew . 

6 Abbreviations used in the article: 

AHR - Antaramian, Hall and Rasin; 

CKM - Cabibbo, Kobayashi and Maskawa; 

CSY - Cheng, Slier and Yuan; 

FCNC - flavor-changing neutral current; 

LHS - left-hand side; 

MSM - minimal Standard Model; 

RHS - right-hand side; 

RGE - renormalization group equation; 

SM - Standard Model; 

VEV - vacuum expectation value; 

1PI - one-particle-irreducible; 

1PR - one-particle-reducible; 

2HDM - two-Higgs-doublet (Standard) Model; 

2HDM(III) - general two-Higgs-doublet (Standard) Model - “type III”. 
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